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ABSTRACT 


* 


Consideration  is  giver  to  a  mathematical  analysis  of  a  special 
case  of  seismic  wave  phenomena.  In  particular,  the  mixed  boundary  value 
problem  of  Love  wave  propagation  in  a  solid  layer  over  a  solid  half¬ 
space  is  investigated  where  the  layer  undergoes  an  abrupt  change  in 
thickness.  Both  the  layer  and  half -space  are  considered  to  be  homogene¬ 
ous  elastic  media.  Theoretical  background  for  the  physics  of  the  problem 
is  provided  by  presenting  the  fundamentals  of  elastic  wave  propagation 
with  specialization  to  the  love  wave  case  and  by  a  statement  of  the 
physical  nature  and  mathematical  form  of  applicable  boundary  conditions. 
Interest  is  focused  on  the  amplitudes  of  the  transmitted  and  reflected 
Love  waves  relative  to  the  magnitude  of  the  excitation,  i.e.,  the  trans¬ 
mission  and  reflection  coefficients,  for  a  range  of  both  the  "strength" 
of  the  discontinuity  (magnitude  of  change  in  thickness  relative  to  layer 
thickness)  and  the  layer  thickness.  The  analysis  employs  scattered 
fields  in  the  form  of  integrals  in  the  complex  plane;  boundary  conditions 
are  applied  to  the  total  fields.  To  satisfy  the  resulting  equations  for 
the  boundary  conditions,  the  nature  of  some  of  the  unknown  coefficients 
in  the  scattered  field  integrals  is  postulated  and  a  function- theoretic 
argument  is  employed  to  determine  these  coefficients.  The  transmission 
and  reflection  coefficients  are  then  extracted  by  a  standard  appeal  to 
the  calculus  of  residues  and  the  energy  contained  in  the  diffracted 
fields  is  evaluated.  Both  displacement  amplitude  and  energy  coefficients 
are  displayed  graphically  as  a  function  of  the  parameters  h  and  H.  The 
results  presented  show  that  the  reflected  energy  constitutes  less  than 

one  percent  of  the  incident  energy  for  all  crust  thicknesses  considered 
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here.  The  amplitude  transmission  coefficient  relating  the  relative 
magnitude  of  the  displacements  of  the  transmitted  and  incident  waves 
is  shown  to  take  on  values  greater  than  unity  for  low  frequency  and 
an  energy  analysis  shows  that  this  behavior  does  not  violate  the  prin¬ 
ciple  of  energy  conservation.  A  comparison  reveals  that  the  results 
shown  here  are  in  close  agreement  with  those  given  by  Knopoff  and 
Hudson  (reference  8).  However,  a  discrepancy  does  exist  for  inter¬ 
mediate  values  of  crust  thickness  and  is  discussed  in  some  detail. 
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INTRODUCTION 

Recently,  considerable  progress  has  been  made  in  the  theoreti¬ 
cal  analysis  of  wave  characteristics  of  various  complicated  structures 
that  approximate  seismic  discontinuities.  Deyond  the  results  that  they 
can  supply,  these  analyses  are  important  both  as  an  indication  of  their 
usefulness  in  examining  other  problems,  and  as  groundwork  for  an  in¬ 
sight  into  more  physically  realistic  situations.  The  problem  considered 
here  is  that  of  Love  wave  propagation  in  a  surface  layer  over  a  half¬ 
plane  substratum.  The  object  of  the  analysis  is  the  determination  of 
the  amplitude  coefficients  and  energy  associated  with  the  transmitted 
and  reflected  Love  waves  which  occur  when  a  Love  wave  is  incident  upon  a 
discontinuity  in  the  layer.  The  topic  is  of  considerable  geophysical  in¬ 
terest  because  it  is  an  approximate  model  of  seismic  wave  interaction 
with  the  continental  margin  and  other  discontinuities  in  the  earth’s 
crust.  Although  the  problem  has  been  considered  by  other  authors,  e.g., 
references  8  and  13,  in  most  cases  these  analyses  lack  graphical  results 
or  include  approximations  which  limit  the  range  of  accuracy. 

A  geometrical  model  for  the  problem  is  shown  in  Figure  1;  a 
fundamental  mode  Love  wave,  propagating  toward  the  right  in  the  surface 
layer,  or  crust,  of  thickness  h,  encounters  a  discontinuity  where  the 
surface  layer  thickness  becomes  h  +  H,  Rigidities,  densities,  and  ve¬ 
locities  of  shear  waves  are  taken  as  |ip  p^,  3^  in  the  surface  layer  and 
*■*2’  ^2’  ®2  in  t*ie  semi_in£ inite  substratum  or  mantle.  The  region  above 
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the  crust  is  taker  to  be  free  space.  It  will  be  shown  later  that,  in 
order  for  a  Love  wave  to  exist  in  this  crust-mantle  structure,  32  and 
must  be  greater  than  ^  and  respectively.  The  origin  of  coordinates 
is  taken  at  the  lower  edge  of  the  step  change  in  crust  thickness,  the 
positive  x  direction  being  the  direction  of  propagation  and  the  z-axis 
pointing  into  the  mantle. 

Unfortunately,  the  geometry  of  Figure  1  leads  to  a  difficult 
boundary  value  oroblenu  In  order  to  facilitate  a  solution,  a  modifica¬ 
tion  of  the  geometry  identical  to  that  used  by  Kane  (reference  7)  is 
applied  to  the  structure.  It  is  felt  that  this  modification  alters  the 
problem  only  slightly  for  small  changes  in  layer  thickness.  The  free 
surface  boundary  x-o,-Hs:zsOis  extended  back  to  x  -  -  ®  as  shown 
in  Figure  2.  In  crder  to  preserve  the  free  surface  of  the  layer  for 
x  <  0,  this  extension  is  performed  in  such  a  fashion  that  there  is  in¬ 
finitesimal  separation  between  the  newly  introduced  waveguide  and  the 
surface  layer  ot  thickness  h,  i.e.,  a  thin  fissure  along  the  negative 
x-axis  that  prevents  transfer  of  radiation  across  that  boundary. 

Since  the  wavelengths  of  seismic  vibrations  associated  with 
the  problem  are  very  long,  it  is  felt  that  the  addition  of  the  new 
acoustic  duct  introduces  at  most  only  minor  perturbations  into  the  pro¬ 
blem.  Indeed  it  is  to  be  expected  that  very  little  energy  would  propa¬ 
gate  into  this  duct  because  the  frequencies  are  in  many  cases  below  the 

cut-off,  and  because  the  duct  is  in  the  backscattered  direction  of  the 
assumed  excitation. 

Finally,  the  media  of  both  crust  and  mantle  are  taken  to  be 
isotropic  elastic  solids.  In  this  case  the  number  of  elastic  constants 
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(which  in  the  generalized  form  of  Hooke's  law  is  36)  degenerates  to  two 
(page  5  of  reference  5),  greatly  simplifying  the  stress-strain  relations. 
In  Chapter  III,  the  case  of  slight  damping  is  assumed,  but  this  is  solely 
for  temporary  mathematical  convenience  and  should  not  be  considered  an 
indication  of  imperfec.  elasticity. 
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II 

FUNDAMENTALS  OF  ELASTIC  WAVE  PROPAGATION 
The  problem  being  considered  concerns  the  propagation  of 
elastic  waves  in  isotropic  media.  Therefore,  an  outline  of  the  theory 
of  motion  in  elastic  solids  and  the  concomitant  equations  of  motion  is 
pertinent.  A  classic  treatment  of  the  material  that  follows  is  found 
in  reference  14. 

Consider  a  point  P  in  a  deformable  body  with  retangular 
coordinates  x,  y,  z  being  displaced  to  a  new  position  x+u,  y  +  v,  z+w. 
It  is  assumed  that  the  displacements  are  small  enough  so  that  second  and 
higher  order  terms  occurring  in  the  stress  and  strain  components  may  be 
neglected.  Then  the  general  form  of  the  stress-strain  component  relations 
becomes  linear  (reference  14). 

From  the  theory  of  elasticity,  the  arrays 


e 
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e 

xy 

e 

xz 
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e 
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yx 

yy 

yz 

e 

e 

e 
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zy 

zz 
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strain  and 

stress  tensors  at  P,  respectively.  Since 

e  ■  e 

t  •  •  • 
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p  ,  ...,  these  arrays  are  symmetrical. 
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If  the  coordinate  system  coincides  with  the  principal  axes,  the  shear 
components  of  stress  and  strain  vanish.  Then  the  deformation  at  P  is 
completely  specified  by  the  corresponding  extensions  (neglecting  higher 
order  t'  rr  s' 


du  dv  dw 

(3)  e-r—  ,  e  -  t-  ,  e  -  r- 

XX  dx  yy  dy  22  02 

and  the  stress  at  P  is  specified  completely  by  the  principal  stresses 

p  ,  p  ,  p  corresponding  to  these  axes, 
xx  yy  z  2 

Consider  a  minute  portion  of  the  matter  enclosing  P  which 
undergoes  deformation.  The  cubical  dialation  6  is  defined  (reference  3) 
as  the  limit,  as  the  surface  area  approaches  zero,  of  the  proportion¬ 
ate  increase  in  the  volume  of  this  matter  and  is  equal  to  the  sum  of 
the  principal  extensions  associated  with  this  deformation: 


(4) 


0  -  e  +  e  +  e 
xx  yy  zz 


du  t  dy  t  dy 
dx  dy  d2 


For  a  derivation  of  the  equations  of  motion,  consideration 

is  given  to  the  stress  components  (2)  across  the  surfaces  of  a  volume 

element  with  uimensions  Ax,  Ay,  and  Az.  The  equations  are  obtained  by 

adding  the  forces  acting  on  the  element  and  the  inertia  terms 
d2u 

-p  — r  AxAyAz,...,  for  each  component.  This  summation  yields: 
dt 


(5) 


.2  dp  dp  dp 

d  u  *xx  ,  yyx  .  *zx 

p  2  "  *  - 

dt  dx  dy  dz 
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(6) 


d2v 

P  ~ 

dt 


dp  dp  dp 

+  J.XL  4.  JLEL 


dx 


dy 


dz 


(7) 


,2  dp  dp  dp 

p  -  w  _  xz  +  vz  +  rzz 

dt2  dx  dy  dz 


In  these  expressions,  body  forces  are  assumed  to  be  absent  and  p  is  the 
density  of  the  medium. 

Assuming  that  the  media  involved  can  be  represented  by  iso¬ 
tropic  elastic  solids,  the  stress -strain  relations  may  be  written  in 
terms  of  two  elastic  constants  (in  this  case  the  Lame'  constants  X  and  p. 
are  arbitrarily  chosen)  in  the  following  manner  (refe^' "*e  b)  : 


(8) 


xx 


yy 


zz 


xe  +  2| i 

xe  +  2p. 


du 

•  p( 

du 

dx 

*Xy 

dy 

dv 

-  Ift( 

dv 

dy 

PyZ 

dz 

dw 

-  n( 

dw 

dz 

Pzx 

dx 

dv 

dx 

dw 

dy 

du 

dz 


Using  (8)  in  *5)-(7)  results  in  the  equations  of  motion  in  terms  of  the 
deplacements  u,  v,  and  w  of  a  point  in  an  elastic  solid: 


(9) 


(10) 


d2i  t  \  4.  .  _2 

p  — 5  -  (X  +  m-)—  +  u 


dt 


d2v  .  Nde  ,  _2 
P  ^2  "  U  +  +  ^  v 


(11)  p  -  (X  +  p)||  +  pV2w 

dt 


2  2 

Note  that  has  been  replaced  by  — r  in  the  left-hand  side  of  (9)-(ll). 
dt*  dt 
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The  justification  for  this  has  already  been  mentioned,  i„e.,  that  the 
second  powers  and  products  which  constitute  .he  difference  between  the 
expressions  are  assumed  to  be  small  -  ignoring  these  products  linear¬ 
izes  the  differential  equations.  The  theory  developed  on  this  basis 
has  been  called  "infinitesimal  strain  theory"  (reference  3). 

It  is  convenient  to  express  the  displacements  in  terms  of  a 
scalar  potential  0  and  a  vector  potential  aS  follows: 

(12)  s(u,v,w)  -  grad  0  +  curl 

Substitution  of  (12)  into  (4)  yields: 

(13)  9  -  V20 


An  expansion  of  (12)  gives  displacement  components  in  terms  of  0  and 
Substitution  of  the  expression  for  u  from  (12)  and  0  from  (13)  into  (9) 
gives : 


d  ,  ,  .  a  ,  Lia  v  _L 

p^?>  ^(p17>  ■»* 


'  W,  1  “ 


.  .i  ,  8  -2.,  .  8  .,2.  8  2, 

-  tt  ♦  V  «  +  I*  s  v  «  ♦  n  *  *3  -  I*  j;  V  i2 


This  equation  and  the  two  others  that  result  from  the  substitution  of 
(13)  and  the  appropriate  component  from  (12)  into  (10)  and  (11)  are 
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satisfied  if  the  functions  0  and  are  solutions  of 


(15) 


V20 


_1_ 

a2  St2  ’ 


B2  dt 


4 

2  ’ 


i  -  1,2,3 


where 


provided  that  p  is  independent  of  the  coordinates  x,  y,  and  z.  Equa¬ 
tions  (15)  are  the  reduced  wave  equations,  and  fhey  indicate  that  the 
existence  of  two  types  of  waves  with  velocities  a  and  0  is  possible  in 
an  isotropic  elastic  solid. 

BOUNDARY  CONDITIONS 

Since  the  problem  deals  with  bounded  media,  some  special 
conditions  expressing  the  behavior  of  stresses  and  displacements  at 
boundaries  must  be  included.  In  particular,  at  the  free  surface  of  a 
solid,  all  stress  components  must  vanish  and,  assuming  that  solid  media 
are  joined  perfectly  at  the  surface  of  contact,  all  stress  and  dis¬ 
placement  components  are  continuous  across  the  interface  between  two 
media. 


LOVE  WAVES 

Under  the  assumptions  of  small  displacements  and  the  absence 
of  body  forces  the  first  section  of  this  chapter  has  shown  that  two  types 
of  disturbances,  governed  by  the  wave  equations  (15)  and  with  velocities 
given  by  (16),  can  travel  in  c  homogeneous  isotropic  elastic  solid. 
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In  addition,  the  previous  section  states  the  conditions  that  the 
stress  components  (8)  and/or  the  displacements  are  subject  to  at 
the  boundaries  ox  the  media. 

The  waveform  of  these  disturbances  is  influenced  by  the 
particular  geometry  of  the  physical  environment.  For  the  stratified 
structure  of  this  problem,  Love  (reference  10)  has  shown  that  one  of 
the  possible  disturbances  consists  of  horizontally  polarized  shear  (SH) 
waves  which,  in  the  notation  of  Figure  1,  means  that  displacements  are 
in  the  y-direction  (v  component).  Although  other  types  of  disturbances 
are  possible  in  this  structure,  the  concern  here  is  with  these  SH  waves 
which  are  called  Love  waves  after  the  person  who  provided  their  explana¬ 
tion. 

-  i<ut 

For  monochromatic  vibrations,  a  time  factor  of  e  can  be 
suppressed.  With  u  *  w  -  o  and  the  displacements  independent  of  the 
coordinate  y,  the  equations  of  motion  (5) -(7)  reduce  to 

(17)  (V2+  k2  )v.  -  0 

81  1 

for  the  layer  and  to 

(18)  (V2+  k2  )v  -  0 

s2 

for  the  substratum,  where 

k  -  7-  ,  kQ  - 

81  sx  &2  e2 


(19) 
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the  shear  wave  velocities  0^  and  02  being  given  by  the  second  expression 
in  (16)  with  p  and  p  subscripted  appropriately.  For  a  crust  thickness 
t  overlying  a  halfspace  and  the  coordinate  system  oriented  as  in 
Figure  1,  with  origin  at  the  free  surface,  assume  displacements  of 
the  form 


(20)  v  -  V  A(y)  cos  (s  z)  eiV3c,  0  a  z  £  t 

l  o,  l  ’ 

(21)  v  z  >  c 

l  o 


vhere 


(22)  »!  -  »Yl-»,/l^')2- 1  -.M-..2 


c  -  phase  velocity  of  the  Love  wave 

and  Vq  is  an  arbitrary  amplitude  constant.  The  boundary  conditions 

require  that  the  stress  p  vanish  at  the  free  surface  z  »  o  and  that 

zy 

Vi  '  <pzy)2 ’  vi  ■  v2  at  the  interface  z  -  t.  By  (8) ,  with  w  -  o, 
(20),  and  (21),  application  of  these  conditions  results  in  A(v)-sec(Sjt) 
and,  for  a  non-trivial  solution,  the  period  equation 


(24) 


tan(s  t)  -  -  -2-* 

1  M1Vl 


The  real  values  of  v  that  satisfy  this  period  equation  are  the 
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U) 

propagation  constants  y  -  —  of  the  Love  wave  (Note  that  -y  are  also 

n  c  *i 

n 

roots  of  (24)  by  virtue  of  its  evenness  in  y) .  Figure  3  shows  a  plot  of 
the  curves  representing  the  right  and  left  hand  sides  of  (24)  as  a  func¬ 
tion  of  the  parameter  s^t.  The  intersections  shown  by  the  circles  de¬ 
fine  the  values  of  y  which  are  real  roots  of  the  period  equation  and 

n 

correspond  to  the  various  possible  modes  in  the  crust,  where  c^  is  the 

phase  velocity  of  the  n—  mode.  In  practical  seismology,  where  much  of 

the  total  energy  is  associated  with  the  smaller  values  of  y^donger 

wavelengths),  the  fundamental  mode  (s^t<tr/2)  is  of  prime  importance. 

For  this  problem,  the  excitation  will  consist  of  only  this  mode.  Although 

the  constructed  solution  will  be  in  general  a  superposition  of  all  possible 

diffracted  modes,  transmission  and  reflection  coefficients  for  the  higher 

order  modes  will  not  be  calculated. 

The  appearance  of  the  radicals  s^^  and  s^  in  the  displacements 

(20),  (21),  in  the  light  of  (22),  (23),  and  (24)  requires  that  k  <v  <k 

“2  n 

or  <  c  <  3^.  Also,  it  can  be  shown  that  0£  <  3^  produces  no  relevant 
solutions.  From  this  fact  and  the  expression  for  0  in  (16)  it  is  estab¬ 
lished  that  the  existance  of  a  Love  wave  is  contingent  upon  ^  ^1  ,,al 

respectively. 

The  phase  velocity  c  is  a  function  of  the  layer  thickness  t. 
Figure  4  shows  this  functional  dependence  by  a  semi- log  plot  of  c/0^  vs. 
the  parameter  yt  for  the  case  ^2^1  "  1.8,  $2^1  *  **29,  where  c  is  the 
first  mode  phase  velocity.  The  parameter  yt  is  a  measure  of  crust  thick¬ 
ness.  As  t  -*  «,  c/8^  -•  1,  i.e.,  as  the  crust  becomes  a  half  plane  the  v^ 
portion  of  the  Love  wave  approaches  a  plane  wave  traveling  with  velocity  0^. 
On  the  other  hand,  as  yt  -*  o,  c/fl^  -  anc*  the  v2  Portion  °f  the  wave 

approaches  a  plane  wave  with  velocity  02« 
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III 

CONSTRUCTION  OF  THE  SOLUTION 

Now  that  the  form  of  the  elastic  disturbance,  the  stress  rela¬ 
tions,  and  the  boundary  conditions  for  the  problem  at  hand  have  been  es¬ 
tablished  it  is  possible  to  assume  a  fundamental  Love  wave  excitation  and 
proceed  with  the  solution. 

For  the  structure  of  Figure  2,  this  excitation  takes  the  form 
(omitting  e"iU,t) 


(25) 


vil  "  Ai(oL>  cos  Cs1(o)*]  e  ;  0  s  z  s;  h 

v  «  e-s2(o)(z-h)  +  icrt 

vi2  e  L  ,  z  >  h 


-  OO  <  X  <  08 


where 


(26) 


A^v)  -  sects  (v)h] 


The  root  «  ■  a  of  the  period  equation 


(27) 


tan£s^(v)h] 


^2  S2 


*1  S1 


(v) 

Tv) 


for  x  c  0  is  the  propagation  constant  tor  this  first  mode  incident  wave. 
As  a  matter  of  referent  this  incident  field  is  assumed  to  have  unit 
displacement  magnitude  at  the  interface  z  ■  h.  Transmitted  and  reflected 
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wave  amplitudes  will  be  compared  with  this  unit  reference. 

It  should  be  noted  tb*t  the  period  equation  (27)  for  x  <  0 
arises  from  the  stipulation  that  solutions  of  the  boundary  condition 
equations  are  non-trivial  and  is  thus  an  essential  aspect  of  Love  wave 
formulation.  Therefore,  the  transmitted  Love  wave  has  a  propagation 
constant  K  -  cu/c’  which  is  a  root  of  the  period  equation  for  the  region 
x  >  0  where  this  equation  is  identical  to  (27)  with  h  replaced  by  h  +  H. 

As  the  initial  step  of  the  procedure  one  seeks  to  add  scatter¬ 
ed  fields  v  . ,  v  .  in  the  crust  and  v  „  in  the  mantle  such  that  the 

si  s3  s2 

total  fields 

(28)  v3  -  vg3,  -  H  £  x  £  0,  all  x 

(29)  Vj  -  vu  *  vsl>  0  <  z  £  h,  all  X 

(30)  v ^  -  vi2  +  vg2,  z  >  h,  all  x 

satisfy  the  homogeneous  wave  equations 

(31)  (V2  +  k2  )  v»  -  0,  -  H  >  z  <;  0,  all  x 

P1  J 

(32)  (V2  +  k2  )  v  -  0,  0  <  z  s  h,  all  x 

h  1 

(33)  (V2  +  k2  )  v,  -  0,  z  >  h,  all  x 

P2  L 

and  comply  with  the  dem.-nds  of  the  boundary  conditions 

<34)  (pzy) i  “  (pzy) 2  f  z  "  h*  a11  x 
(35)  v^  -  v2,  z  »  h,  all  x 
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(36) 

(p  ),  ■ 

rzy  3 

0, 

z  - 

-  H 

,  all  x 

(37) 

■ 

0, 

Z  " 

0, 

x  <  0 

(38) 

(V1  ■ 

0, 

z  - 

0, 

x  <  0 

(39) 

(p  ),  - 

rzy  3 

(p  ) 

rzy 

1» 

Z  " 

0,  x  >  0 

(40) 

V  «■  V 

'3  r 

Z  " 

0, 

X  > 

0 

Also,  except  for  the  incident  fields,  v^,  v^,  and  v^  should  represent  out¬ 
going  waves  and  be  bounded  and  continuous  at  the  origin  or  edge  of  the 
fissure. 

Arbitrary  solutions  of  the  reduced  wave  equations  (31)-(33)  can 
be  written  as  follows: 


(4l)  VS1  ■  1ST  I  [ 


C(v)  sin  [s^(z-h)]  d(v)  cos  [s^z-h)] 


s,  cos  s.h 
c  1  1 


s^  sin  Sjh 


du 


,.n\  1  f  F(y)  -s.(z-h)  +  iyx, 

(42)  vs2-toTj  6U)e  2  iV 


....  IP  B(0)  C0S  4„ 

(43)  vs3"  *  5iJ  - sTTiiTsTH -  e  d» 

C  1  1 


where 


(44)  6(v)  -  s1  sin  Sjh  -  —  s2  cos  Sjh 


and  s^,  s^  are  given  by  (26).  The  coefficients  B(y) ,  C(y),  D(y)  ,  and  F(y) 
are  unknown  functions  of  v  and  C  is  some  contour  in  the  complex  y-plane. 
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At  first  it  seems  that  the  form  of  (41) -(43)  introduces  undue  complexity 
into  the  mathematics  but  it  will  be  seen  that  these  expressions  lead  to 
simple  boundary  condition  equations,  thus  facilitating  solutions  for  the 
unknown  coefficients. 

It  is  mathematically  convenient  to  assume  a  small  imaginary 
i  € 

part,  proportional  *'o  e  ,  in  the  propagati  a  constants  k  ,  k  and  in 

B1  e2 

the  roots  a,  K  associated  with  the  regions  x  <  0  and  x  >  0  respectively; 
that  is , 


(45) 


e 


ie 


o-|a|ele  ,  K-|K|ele 


This  assumption  causes  the  fields  to  possess  a  slight  attenuation  in  the 

direction  of  propagation;  in  the  final  analysis  e  will  be  set  equal  to  zero. 

Vertical  branch  cuts  from  the  branch  points  ±  k  are  so  chosen  that  the 

P2 

cut  v-plane  shown  in  Figure  5  represents  that  sheet  of  the  Riemann  surface 
for  which  Re  /  y2-  k2  >  0. 

e2  - 

By  (8)  with  p^  ■  p^z,  application  of  the  boundary  conditions  (34)- 
(40)  to  the  total  fields  (28)-(30)  yields,  from  (34); 


(46) 


2M  J[ 


^(v) 

cos  s^h 


^2s2 

6(v) 


] 


iyx 


du  ■  0,  -00  <  X  <  « 
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from  (35) ! 


(47) 


-L-  f  f  ,  .  liul  1  eiuX  du  -  o 

2tri  .)  l  s.sins.h  6(u)  Je 
cl  i 


,  -OO  <  00  < 


The  form  of  is  such  that  (36)  is  identically  satisfied. 
From  (37) : 


(48) 


^  B(u)  eiuX  du  -  0,  x  <  0 


from  (38) : 


(49)  1ST  J  ^i  [c(u)  +  D(u)]  *iUX  du  "  °*  x  <  0 


from  (39): 


(50) 


2wT  J  "  ^(v)  “  ®iUX  du  -  0,  x  >  0 


Finally,  the  condition  that  expresses  continuity  of  displacement  at 
z  -  0,  x  >  0  includes  a  non-zero  contribution  from  the  incident  Love 
vc.ve.  For  x  positive,  this  excitation  can  be  written 


(50  *11- ITT  J  Al(a)cos[sl(a)l]  J-x 


u-a 


du 


(52) 


1  P  i„x  . 

Vi2-2«J  - ^ -  e  J» 

c 
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by  the  C  hy  Integral  Theorem*  The  boundary  condition  (^0)  yields: 


,  .  r  B(y)cot  s.H 

«»^r.  II- — : — L 


2tri  J  L 
c 


A^(a)  C(y)tan  s^h  -  D(y)cot  s^-  4  x 

-  + - -  e  u  dy  “  0, 

v-  a  s4  J 


for  x  >  0. 


A  beginning  toward  satisfying  the  above  integral  equations 

(46) -(50)  and  (53),  is  found  by  postulating  the  nature  and  behavior  of 

some  of  the  coefficients  in  certain  regions  of  the  y-plane.  Equation  (48) 

is  satisfied  if  B(y)  is  analytic  for  Imy  <  Jk  |  sin  e  and  is  of  the  order 
-cl  .  ,.c  ..  02  . 

v  (tj>  0)  as  |u|  -*  »  in  the  specified  half-plane  of  analyticity.  Then 

the  integral  will  vanish  by  Jordan's  Lemma  (reference  9).  This  postu¬ 
lated  behavior  is  denoted  by  writing  B(y)  as  B~(y).  Equations  (49)  and 
(50)  are  satisfied  if 


C(u)  +  D(v)  -  B  (v) 


Consider  next  (46)  and  (47)  which  hold  for  all  x.  Equating  the  integrands 
to  zero  and  solving  for  C(y)  and  D(u)  yields: 

4,  s,cos  s  h 

(55)  C(”>  ■  -  sj  Hv) 

s  sin  s  h 

(56)  D(v)  -  F(w) 

Inserting  the  expressions  into  (54)  with  6(y)  given  by  (44)  results  in: 

Kv)  -  r‘(u) 


(57) 
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nus 


,  F(u)  has  the  same  behavior  as  B  (y)  for  Imy  <  jk  |  sin  e, 

°2 

i.e.,  F(y)  A  F  (y) ,  The  last  equation  to  be  satisfied  is  (53). 

With  F(y)  -  B  (y)  and  C(y)  and  D(y)  given  by  (55)  and  (56),  (53)  then 
becomes 


(58) 


S7  I  G(l,) 


iyX, 

e  dy 


0. 


x  >  0 


where 


(59) 


A. (a)  r  cot  s  H  |i  s0  sin  s  h  cos  s.h 

G(rf-  -  -tr *  [-  — r1-  -  --5  -  -st4-  »-<v> 


1 


M1S1  6(y) 


6(y) 


This  last  boundary  condition  equation  is  satisfied  if  G(y)  is  analytic 

:2, 

I  sin  g  ana  is  nr  rnp  nrop.r  11 
0  2 


-e, 

for  Imy  >  -  M  sin  e  and  is  of  the  order  y  ie  >0)  as  Jyj  -*  ®  in 


this  half  plane.  This  behavior  is  denoted  by  writing  G(y)  as  G  (u). 

Now  the  problem  becomes  that  of  solving  for  the  unknown  co¬ 
efficients  of  the  scattered  fields,  particularly  B  (y)  which  has  the 
postulated  behavior. 

Since  (59)  is  the  only  expression  available  for  the  extrac¬ 
tion  of  the  unknowns  it  becomes  necessary  to  introduce  a  function- 
theoretic  argument  based  on  the  Wiener-Hopf  technique  (reference  12). 
For  convenience  let 


(60) 


£(y) 


cot  SjH  •“l2S2  sin  si^  cos  s]k 
s^  l^s^  6(y)  6(y) 


Provided  that  it  has  certain  properties,  £(v)  may  be  expressed  as 
(references  2,  11,  12,  and  16,) 
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(61) 

£(v)  -  2— 

Axl 

CT 

(v) 

where 

1. 

CT  +  (v) 

is 

regular  and  zeroless  for  Imv  > 

-  Ikq  1  sin  e* 
92 

2. 

CT  (v) 

is 

regular  and  zeroless  for  Imv  < 

|k  |  sin  e. 

d2 

3. 

k  +  (v)i,| 

k 

(y)|  lie  between  positive  bounds  in  their 

respective  domains  of  regularity  as  |y|  becomes  large. 
Consequently,  (59)  can  be  rewritten,  on  dividing  through  by  a  +  (v),  as 


(62) 


Ai(a) 


G_M  +  . 

cr  (v)  (v-a)  a  (v) 


- 

CT  (v) 


Subtracting 


A^q) 


from  each  side  gives 


(63) 


(v-a) 

ct  (a) 

-  Cv) 

At(a)  . 

1  l 

]  . 

B-(„  Ai(a) 

+  (v) 

v-  a[ 

r  (v)  ct  (a 

J 

ct  (v)  (v-a)  CT+(a) 

An  inspection  of  (63)  reveals  that  the  left  hand  side  is 

regular  for  Imv  >  -  |k  |  sin  c  and  the  right  hand  side  is  regular 

for  Imv  <  |k  j  sin  e.  In  addition,  because  of  the  assumed  growths  of 

♦  .2 

G  (v)  and  B  (v)  for  large  Jy|,  these  functions  tend  to  zero  as  |v|  -*  » 
in  their  respective  half  planes. 

By  the  previous  paragraph  and  the  overlapping  regions  of 
regularity  shown  in  Figure  5,  each  side  of  (63)  is  the  analytic  continua¬ 
tion  of  the  other  (reference  15)  and  thus  defines  an  entire  function. 
Invoking  the  Liouville  Theorem  (reference  4)  in  the  light  of  the  prescribed 
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behavior  at  |p|  -  each  side  of  (63)  can  be  equated  to  zero,  giving 


(64) 


(65) 


G  (v) 


B‘(v)  - 


A  (a)  _  +( 

_i -  rg_  (v)  .1 

v  -  a  .  "  AJ 

a  (a) 

A.(a)  ct"(p)  o+(p)  A. (a) 

(v-a)  a  ia)  a+(a)  (p-a)  £(p) 


ia 


The  diffracted  fields  can  now  be  specif '.ed  completely  vi 

substitution  of  (65)  and  (55). (57)  into  (41). (43),  and  the  total  fields 
(28) -(30)  become 


l  r  a  (p)  A  (a)  cos  [s  (z  +  H)] 
(oh)  -  - — T  -  — - - - 1 _ 

J  2ffl  J  +.  .  ,  . 


ca(a)  (p-a)S(p)s1  sin  s^ 
A. (a) 


ipx  . 

e  du,  -H  s  z  £  0 


(67)  vx  vu+  2ni  (y_c^i;(i?)6(v)['cos[si(z-h)]  -  ~~  sin[s1(z-h)]]eivxdv, 


0  s  z  s  h 


(68)  v  -  v  +  JL  I*  2_IvI  _ S  ^ 

2  12  2ni  c  a+( a)  (u-a>S:(v)«(v) 


A. (a)  -s  (z-h)+ipx 

e  dp. 
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IV 

THE  TRANSMISSION  AND  REFLECTION  COEFFICIENTS 
Evaluation  of  Total  Fields 

The  integral  expressions  for  the  total  fields  (66) -(68)  can,  at 
least  in  principle,  be  evaluated  by  a  standard  method.  For  the  Love  waves, 
depending  on  whether  the  transmitted  (x  >  0)  or  reflected  (x  <  0)  field  is 
being  considered,  the  contour  C  along  the  real  axis  is  closed  in  the  appro¬ 
priate  half  plane  by  an  infinite  arc  with  due  regard  being  taken  not  to 
cross  the  branch  cuts  emanating  from  the  branch  points  ±k  .  By  Jordan's 
Lemma  the  integral  along  the  infinite  arc  makes  no  contribution  so  that  by 
the  Residue  Theorem  (reference  4)  the  closed  contour  integral  can  be  ex¬ 
pressed  as  follows: 

(69)  ]  -  +  F  branch  cut  *  ±  2iri  £  residues, 

closed  c  J 
contour 

The  branch  cut  integral  can  be  evaluated  by  the  method  of  saddle  points. 

The  resultant  field  is  the  SH  cylindrical  body  waves  that  are  scattered 
into  the  mantle  and  will  not  be  discussed.  Therefore,  by  (69)  the  scattered 
surface  waves  of  interest  are  obtained  as  the  residue  fields  of  the  closed 
contour  integral. 

For  the  transmitted  field,  application  of  Jordan's  Lemma  re¬ 
quires  that  the  contour  be  closed  in  the  upper  half  plane.  Letting  e 
shrink  to  zero,  the  contour  C  is  deformed  below  singularities  for  Rev  >  0 


22 


and  above  these  singularities  for  Rev  <  0  as  shown  in  Figure  5. 

This  contour  and  the  upper  half-plane  with  both  Reimann  sheets  displayed 
is  shown  in  Figure  6.  The  period  equation  for  the  region  x  >  0  is  in¬ 


cluded  in  E(v)(see  equation  3,  Appendix  A)  and  produces  a  pole  of  the  in- 

I  2  2 

tegrands  of  ( 66) - (68)  at  v  -  ±  K  provided  that  the  radical  s  -  v  -k 

2  V 

is  positive  real.  Therefore,  the  poles  ±  K  lie  on  the  sheet  that  ex- 

I  2  2 

tends  from  the  upper  left  to  lower  right  of  Figure  6  where  Re  v  -k  >  0 

_  V  Bn 


I  2  2 

general,  and  where  /  v  -k  is  for  positive  real  for  v  real 

'  P2 


greater  than  k  . 


The  integrand  in  (66)  has  no  other  poles ,  but  those  in  (67) 


and  (68)  have  poles  at  v  -  ±  a  provided  that  s^  has  the  same  sign  as  before, 
as  can  be  seen  from  the  period  equation  (27)  for  the  thinner  layer.  That 


these  poles  at  ±  a  exist  is  shown  as  follows:  Expressing  6(v)  as 


r^2  s2 

6(v)  -  -  s.  cos  s,h  -  - 

1  1  L^i  si 


tan  s. 


from  (44)  and  considering  the  product  (v-a)E(v) 6(v) ,  one  sees  that  there 

is  a  cancellation  of  the  period  equation  in  h  between  the  denominator  of 

E(v)(see  equation  3,  Appendix  A),  and  6(v) ,  leaving  the  linear  factor  (v-a) , 

thus  giving  a  pole  at  v  ■  a.  For  the  pole  at  v*  -  a,  consider  E(v)  as  the 

quocient  g .  Then, 
a  (v) 


(71) 


+/\  A. (a)  “/  \  A  ( a) 

SLivI  ,  i _ o  Cv)  i 

\  (v-a)E(v)6(v)  +,  x  (v-a)6(v) 

a  (a)  a  (a) 


and  6(v)  contains  the  pole  at  v 


-  a  . 
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The  pole  at  v  -  -  a  is  associated  with  the  reflected  wave. 

In  this  case,  the  application  of  Jordan's  Letr  >0  requires  that  the  contour 
be  closed  in  the  lower  half  plane.  This  suitably  indented  contour, 
both  Riemann  sheets  of  the  lower  half  v-plane,  and  the  pole  at  v  ■  -  a 
are  shown  in  Figure  7.  The  sheets  in  this  figure  are  oriented  as  in 
Figure  6,  i.e.,  Re  s^t  0  on  the  sheet  extending  from  upper  left  to  lower 
right.  This  orientation  does  not  imply  that  the  two  sheets  have  the  en¬ 
tire  branch  cut  in  common.  In  fact,  only  the  branch  point  is  common  to 
both  sheets. 

It  must  be  noted  that  no  mention  has  been  made  of  the  singulari¬ 
ties  of  a+(v)  and  a"(v).  Although  both  these  factors  contain  singularities, 
they  are  used  in  the  integrals  so  that  these  singularities  are  not  included 
within  the  contours  of  Figures  6  and  7.  The  transmitted  field  integrand 
contains  a+(v)  which  is  regular  and  zeroless  in  the  upper  half  v-plane, 
which  is  the  region  enclosed  by  the  contour  of  Figure  6.  For  the  re¬ 
flected  field,  the  integrand  may  be  altered  by  use  of  (71)  and  the  term 
ct  (v)  presents  no  singularities  within  the  contour  of  Figure  7. 

The  specification  of  contours  and  location  of  poles  makes  it 
possible  to  extract  the  residue  fields.  Noting  that  the  residue  associated 
with  the  pole  v  ■  +  a  serves  only  to  cancel  the  incident  Love  wave  in  the 
region  x  >  0  (which  is  to  be  expected  since  the  Love  wave  cannot  propagate 
undisturbed  past  the  discontinuity),  the  transmitted  fields  are 


(72)  vtl-Tl(tcB 


1 


,k  ,h,H){cos[s  (K)(z-h)]- 
°2  1 


u2s2(K) 


sin[s1(K)(z-h)])eiKx, 


0  s:  z  <;  h 
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-r  i  u  u\  -S-(K)(z-h)  +  iKx  w  , 

(73)  vt2“  T2{\*  h»  H)  e  2  ,  z  >  h 


(74) 


vt3  "  Vk0  »  k02»  h»  H)  cos  tsi(K>(z  +  H>3  elKX» 


•H  s  z  s  0 


where  the  transmission  amplitude  coefficients  T^,  T2,  T^,  can  be  calculated 
as 


(75) 


(76) 


T  .T  _ Ai(oJ 

1  2  o' \o)  (K-o)  SC»)  I  600 


V"K 
Aj(a) 


T  -  -  ?  W  - 

J  a+( a)  (K-a)  4-  E(v)  |  s^.K)  sin  [«1(K)H] 

u*K 


with  Sl(K)  -  Jk*  -  K2  ,  s2(K)-^2  -  k2 


The  reflected  fields  are 


n2s2C-J 


(77)  v  R. (k  ,  k  h,  H) {cos  [s. (-a) (z-h)3-  ; — r  sin[s. (-a)(z-h)]}e“iaX, 

rl  1  8i  @2  1  UjS^i-a)  1 


0  is.  z  <.  h 


(78)  vr,  -  R,(k  ,  kQ  ,  h,  H)  e-s2(“a3(z”h-)  -°oc>  z  >  h 


where,  by  virtue  of  equation  15,  Appendix  A,  the  reflection  amplitude 
coefficients  R^,  R2  are! 


(79) 


*1  -  R2  ■  2 


Jz 


-a-ijk  Jsin  e 
P 


Ai(a) 


a+1  |ke2|sin  6  [oW  2a^6(v)| 


w-a 
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with  s^(-a)  m  Jk 2  -a2  ,  s2(-a)  -  Ja'-  k2  . 


The  derivatives  in  (75),  (76),  and  (79)  are  as  follows! 


(80)  £  S(„)  .  ±  {.[ 


2  h 

,,  2  „  u  vh  —  cos  s.h-sin  s.  — 

VH  COS  S.H  COt  S.H  _  |i  _  s.  1  1 

=  +  ~  ] +  IT  [  “  “  ~ + 


Y, 


Sj  6(v)/v 


2 


S2  sln  slh  C*1  £  Hv)-  ^  6<”>:l, 

si  y2C6(o)32/u2 


-r 


vh  sin  s^h 
Hv)/v 


6(v)  cos  s,h 


[6(v)]2/v2 


*=]} 


d  f  1  ^oS2 

(81)  —  6(v)  -  -  4  yh  cos  s.h  +  —  sin  s.h  +  -  vh  sin  s.h  + 

dv  1  1  1  li^s^  1 

(i_  cos  s.h  , 

♦  -2  - i-  L 

“l  ^2  J 


where  y^  y2  are  given  in  (22)  and  (23). 


Energy  Considerations 

For  a  Love  wave  in  the  fundamental  mode,  the  rate  of  energy 
which  passes  a  plane  perpendicular  to  the  x-.xxis,  per  unit  area,  per  unit 
time  is  (reference  13) 


(82) 


E  «±n 


dv 

3x 


dv 

dt 


where  the  minus  sign  is  used  for  propagation  in  positive  x  direction  and 
the  plus  sign  for  propagation  in  the  opposite  direction.  For  a  certain 
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with 


s^C-a)  -  Jk^  -a2  ,  s2(-o)  -  J<x  -  k2  . 


The  derivatives  in  (75),  (76),  and  (79)  are  as  follows; 


(80)  i  S<„)  -  \  {-[ 


vH  cos  s.H  cot  s,H 

- - — L_  +  — 


*2  .  h 

vh  —  cos  s,h-sin  s.  — 

1  Vo 


iH  -I  /2  r  9i  1 
J  L 


s^  6(v)/v 


s2  sin  sih  [VX  ~  ecu)] 

si  Yj[6<»)]2/#2 


]-[ 


vh  sin  Sjh 
6(v)/v 


^  fi(v)  cos  s , h 


[6(v)]2/v2 


L  } 


(81)  6(v)  -  -  |  vh  cos  s jh  +  —  sin  s^h 


^2S2 

+  — vh  sin  s.h  ♦ 
¥i  1 


\i~  cos  s.h  , 
* - 1-\ 


‘i  v2 


where  y  ,  y2  are  given  in  (22)  and  (23). 

Energy  Considerations 

For  a  Love  wave  in  the  fundamental  mode,  the  rate  of  energy 
which  passes  a  plane  perpendicular  to  the  x-axis,  oer  unit  area,  per  unit 
time  is  (reference  13) 

(82)  E  “±(i  ~  |^ 

ax  at 

where  the  minus  sign  is  used  for  propagation  in  positive  x  direction  and 
the  plus  sign  for  propagation  in  the  opposite  direction.  For  a  certain 
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period  T-2it/u>  the  energy  flow  is 

(83)  E  -  f  Edt 

J0 

and  the  total  flux  across  a  vertical  plane  of  unit  breadth  in  a  period  T  is 

<**- 

(84)  [T.F.]  -  j  E  dz 

*0 

Expressing  the  incident  field  as  the  real  part  of  (25)  with 
the  time  factor  e"iUlt  included  yields: 

(85)  v  .  -  A  cos  s.z  cos(-u>t  +  ax) ,  0  s  z  S  h 

ii  i  i 

(86)  vj2  -  e’S2^z“h^  cos(-u)t  +  ox) ,  z  >  h 

where  and  are  understood  to  be  s^a)  and  s2(a),  respectively,  and 
Aj  -  A^a)  (see  (26)).  The  energy  rate  per  unit  time  and  area  of  the 
incident  wave  is  then 

2  2  2 

(87)  E^  -  A^  cos  s1  z  sin  (-u>t  +  ox) 
for  the  crust,  and 

-2s2(z-h)  2 

(88)  Ei2  *  aU>^2  e  sin  (-u>  t  ♦  ax) 


for  the  mantle. 


Integrating  over  a  period  T  gives 


!  "Ilf  t  ftfi&fitiltf  ii'HI  H  i  -‘till  *('•  iftitHltl**1  * 
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Finally,  integration  over  the  vertical  coordinate  z  results 
in  an  expression  for  the  total  incident  flux  crossing  a  vertical  plane 
of  unit  breadth  in  time  T: 

0  .  sin  2s  h 

(91)  [T.F.]n-H1.*Ja[|*— <-*-] 

(92)  [X.F.],2  -  J  ^ 

(93)  Cl.F.3inc.  -  "  tT.F.Di2 


The  expression  (77)  for  the  reflected  field  in  the  crust  can 

be  simplified  by  substituting  the  right  hand  side  of  (27)  with  v  -  -  a 

P2s2(-a) 

for  the  quantity  — -  ^  ^  .  This  substitution  simplifies  (77)  to 


vrl  '  Rl\COS  V  e 


-iax 


A  comparison  of  (94)  and  (78)  with  (85),  respectively,  with  real  parts 
considered,  -  R^,  s^-a)  -  s ^(a),  and  s2(-a)  ■  S2^»  shows  that 

the  reflected  fields  differ  from  the  incident  fields  by  the  coefficient 
Rj^  2  and  the  sign  of  a.  The  sign  is  accounted  for  by  the  sign  selection 
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in  (82).  Therefore,  it  is  seen  that  the  reflected  energy  is  equal  to  the 

2 

incident  energy  times,  say,  »  or 

(X  F.] 

(95)  - SSLUm  r2 

[T.  F.].  2 

me. 

This  concise  expression  (95)  is  true  only  for  the  case  where 
the  incident  and  reflected  fields  in  question  are  the  same  mode. 

The  analysis  (82)-(84)  as  applied  to  the  transmitted  fields 
is  facilitated  by  agaki  employing  the  period  equation  for  the  region 
x  >  0  to  simplify  (72)  to 

(96)  vu  -  TjAjCos  [Sl(z  +  H)]  eiKx 
where 


(97)  At  -  At(K)  -  sec  [e1(K)(h  +  H)] 


The  end  result  of  this  energy  cnalvsi-  on  (96),  (73),  and  (74)  is  as 
folic  s : 


(98)  [T.  r.]tl  -  Hjtt  ij  A*  Kf  |  +  i 


sin[2  s^(h  +  H)]-sin  2  s,H] 


1 


1 


(99)  [T.  r.]tl  -  i  V  t*  f. 


(100)  [T.  7.3t3  -  ^  TT  T2  K{  |  +  £ 


sin  2s 
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and  the  total  flux  transmitted  is  equal  to  the  sum  of  (98), (99),  anti 
(100).  In  these  transmitted  flux  expressions,  s^,  are  understood 
to  be  s^(K),  s^(K) ,  respectively. 

Of  interest  here  are  the  energies  contained  in  the  reflected 
and  transmitted  Love  waves  as  a  fraction  of  the  energy  in  the  incident 
Love  wave.  This  ratio  for  the  reflected  wave  is  given  by  (95).  The 
rati ?  for  the  transmitted  wave  is  readily  obtainable  from  the  above  ex¬ 
pressions  and  is  not  written  here  due  to  the  space  involved. 


Graphical  Results 

Up  to  this  point,  the  theory,  its  application,  and  the  re¬ 
sulting  expressions  for  amplitude  coefficients  of  and  flux  contained 
in  the  various  fields  can  have  at  most  a  pseudo- importance  and  conjectured 
validity.  It  may  be  impressive  to  display  the  analysis,  concluding  with 
a  peremptory  statement  such  as  "these  expressions  are  the  — ",  but  the 
work  is  neither  practical  nor  believable  until  the  magnitudes  of  the 
coefficients  and  energy  ratios  are  known  as  a  function  of  the  physical 
parameters  h  and  H. 

The  proceeding  graphs  display  the  quantities  that  are  felt  to 
be  of  prime  importance.  All  major  calculations  were  performed  on  the 
IBM  digital  computer  at  the  University  of  Rhode  Island.  Numerical  com- 


^2 

putation  requires  that  the  values  of  —  and  —  be  specified.  For  the 

^1  S1 

purpose  of  future  comparisons  with  other  authors,  the  values  stated 
in  Figure  4  are  used. 

The  calculations  of  T^ ,  i  -  1,2,3,  i  R. ,  i  -  1,2  requires 
that  consideration  be  given  to  magnitude  only,  phase  being  ommitted,  for 
the  expressions  representing  these  coefficients  contain  terms  which  are 
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is  in  the  shadow  zone  of  the  incident  wave.  Also,  for  the  trans¬ 
mitted  wave  in  the  fundamental  mode^s^CK)  (h+H)  <  tt/2,  so  that  one 
is  inclined  to  believe  that  s^(v)H  will  remain  less  than  tt/2  at  least  for 
fairly  small  H.  This  would  restrict  the  value  of  n  in  (102)  to  zero, 
i.e.,  the  propagation  would  be  analogous  to  the  TEM  electromagnetic  wave 
in  a  parallel  plate  waveguide.  Therefore,  it  is  felt  that  any  propaga¬ 
tion  along  the  duct  for  small  H  will  constitute  only  a  very  s^ail  per¬ 
centage  of  the  energy  associated  with  the  original  problem  of  Figure  1. 

Indeed,  no  unusual  behavior  is  observed  in  T0  for  k  H  >  3,  although 

1  31 

there  is  no  criterion  to  determine  either  the  accuracy  of  the  results 

for  large  discontinuities  or  the  dividing  line  between  valid  and  invalid 

ranges  k  H. 
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Note  that  |T2|  takes  on  values  greater  than  unity  when  the 

crust  thickness  for  x  <  0  is  small  (k  h  -  .05").  However,  this  unusual 

01 

behavior  is  justified  by  a  plot  of  the  -atio  of  transmitted  to  incident 

flux  TFT/TFI  shown  in  Figure  10  for  the  smaller  layer  thicknesses  of 

Figure  8.  This  figure  verifies  the  fact  that  coefficients  larger  than 

unity  are  not  necessarily  in  violation  of  conservation  of  energy.  In 

fact,  Knopoff  and  Hudson  (reference  8)  have  chown  that  the  rate  of  energy 

transmission  in  a  Love  wave  is  inversely  proportional  to  layer  thickness 

for  small  frequencies  (k  h  -  ^  -*  0) .  It  follows  that  a  coefficient 

°1  B1 

associated  with  transmission  from  a  thin  layer  to  a  slightly  thicker  one 
which  is  greater  than  unity  is  consistent  with  energy  conservation  prin¬ 
ciples  for  this  range  of  small  frequencies. 

Figure  10  indicates  that  for  k  h  >  2.040,  the  rate  of  change 

ei 
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of  the  transmitted  to  incident  energy  ratio  is  at  first  rapid  and  then 

diminishes  to  an  almost  constant  value. 

Since  IR^I  is  small,  the  normalized  reflected  energy  is 

negligible  (<  1  percent)  and  will  not  be  displayed. 

Knopoff  and  Hudson  have  presented  an  approximate  solution  to 

the  problem  dealt  with  here  (reference  8).  Their  results  are  plotted  as 

a  function  of  the  crust  thicknes  k  (h+H)  fot  x  >  0  with  the  ratio  of 

1  h 

crust  thickness  for  x  <  0  to  crust  thickness  for  x  >  0  (r—)  held  con- 

u+H 


stant.  The  incident  Love  wave  in  their  paper  has  unit  amplitude  at  the 
free  surface  of  the  crust  and  the  coefficient  calculated  is  the  free 
surface  amplitude  of  the  transmitted  Love  wave.  In  this  thesis  the 
equivalent  quantity  is  readily  obtained  as  Jt^|  divided  byA^cO,  and 
will  be  denoted  by  |tt|.  The  method  of  plotting  is  as  follows:  The 
growth  of  the  crust  is  such  that  the  ratio  of  thin  to  thick  layers  re¬ 
mains  constant,  and  the  curves  are  shown  as  a  function  of  the  crust  thick¬ 
ness  k  (h+H)  for  x  >  0.  Figure  11  is  a  comparison  of  the  results  of  Knopoff 
^1 

and  Hudson  and  this  author  for  ratios  of  incident  to  transmitted  layer 

thickness  of  0.9,  0.8,  and  0.7.  Note  that  the  results  agree  with 

10  percent  for  k  (h+H)  >5,  and  that  the  limiting  values  as  k  (h+H)  -•  0 
01  Bl 

are  coincident.  For  k  (h+H)  >  9,  the  curves  become  asymptotic  to  the 

B1 

dashed  lines  shown  in  the  figure. 

In  order  to  profit  from  this  comparison  it  is  necessary  to  be 
aware  of  any  assumptions  and  the  corresponding  errors  involved  in  the 
method  of  knopoff  and  Hudson.  Only  then  can  the  comparison  become  a  gauge 
to  measure  the  quality  of  the  results  of  this  thesis. 
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First  of  all,  the  displacements  along  the  aperture  plane  x  -  0, 

0  £  z  £  *  were  assumed,  by  Knopoff  and  Hudson,  to  be  thi.se  due  to  the  in¬ 
cident  wave  alone.  This  assumption  is  accurate  as  long  as  the  discontinuity 
is  negligible  with  respect  to  the  wavelength  of  incident  field.  A  small 
or  "weak"  discontinuity  neither  converts  appreciable  energy  to  SH  body  waves 
nor  needs  to  create  evanescent  modes  to  effect  the  transition  from  incident 
to  transmitted  Love  waves.  Also,  from  Figure  9,  the  reflected  wave  is  too 
small  to  introduce  appreciable  error.  Therefore,  the  dashed  line  curves  of 

Figure  11  can  be  considered  very  accurate  as  k  (h+H)-*  0,  for  k  H-*  0 

gl  “l 

simultaneous ly. 

The  quantitative  effect  of  the  above-mentioned  assumption  for 
large  crust  thickness  is  not  so  easily  determined.  Of  all  the  distur¬ 
bances  associated  with  the  problem,  i.e.,  reflected  wave,  SH  body  waves, 
and  evanescent  modes,  only  the  first  can  now  be  considered  negligible.  It 
is  suspected  that  there  is  some  conversion  to  body  waves,  but  the  magni¬ 
tude  and  direction  of  these  body  wave  displacements  are  unknown.  The 
existance  of  evanescent  modes  is  also  suspected,  and  again,  the  magni¬ 
tudes  of  these  modes  are  unknown.  An  indication  of  a  smooth  transition 
from  incident  to  transmitted  Love  wave  is  obtained  from  Figure  4.  This 
figure  shows  that  the  Love  wave  phase  velocity  is  practically  invariant 
vith  crust  thickness  for  large  thicknesses,  inferring  that  the  Love  wave 
is  relatively  unaffected  by  the  discontinuity.  Nevertheless,  the  above 
statements  lead  one  to  be  skeptical  about  the  accuracy  of  Knopoff 's  and 
Hudson's  results  for  k  (h+H)  >6,  On  the  other  hand,  it  is  felt  that 

h 

the  dashed  curves  of  Figure  11  are  a  fair  approximation  to  the  desired 
coefficient  for  large  crust  thickness. 


34 


Finally,  the  assumption  under  consideration  is  probably  least 
valid  for  intermediate  crust  thickness.  In  addition  to  the  existance  of 
evanescent  and  body  wave  disturbances  al!  ig  the  aperture  plane,  Figure  4 
shows  that  the  phase  velocity  changes  considerably  for  a  sm.il  change  in 
crust  thickness,  thus  indicating  a  rough  and  complex  transition  from  in¬ 
cident  to  transmitted  wave.  These  conditions  imply  that  the  intermediate 
crust  thickness  coefficients  of  Knopoff  and  Hudson  are  only  rough  approxi¬ 
mations  . 


The  second  assumption  of  Knopoff  and  Hudson  is  that  the  dis¬ 
placements  along  the  surface  x-0,  -  H  £  z  0  are  negligible.  This  heuristic 
judgement  is  felt  to  be  valid  for  at  least  the  small  discontinuities  asso¬ 
ciated  with  Figure  11. 

The  above  discussion  and  the  comparison  of  Figure  11  establishes 
the  accuracy  of  the  calculated  values  of  this  work  as  k  (h+H)  -  0  and 
provides  a  measure  of  the  accuracy  of  these  calculated  values  for  large 


k  (h+H). 
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Figure  12  shows  the  coefficient  JT-J  in  the  same  manner  and 


for  the  same  values  rrTT  as  Figure  11.  The  limiting  values  as  k  (h+H)-«  0 

h  w  p  ^ 

are  the  same  as  those  for  |tt|  curve  of  Figure  11,  finally  approaching  the 

dashed  lines  shown  for  k  (h+H)  >  9. 

01 


Again,  the  validity  of  the  values  of  |TT |  and  j T^ |  in  Figures  li 

and  12,  at  least  in  terms  of  energy  conservation,  can  be  established  by  a 

plot  of  the  fraction  of  energy  contained  in  the  transmitted  Love  wave  vs. 

k  (h+H).  Figure  13  shows  that  the  coefficients  JttJ  and  |T  |  do  satisfy 
ei  2 

these  energy  limitations  and  that  the  portion  of  energy  transmitted  approaches 

a  "onstant  value  for  large  crust  thickness. 
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A  final  and  interesting  comparison  with  the  results  of  Knopoff 
and  Hudson  is  shown  in  Figure  14.  Calculations  were  performed  on  [ TT [ 

1  ?— assumed  to  have  the  value  unity,  and  the  results 


with  the  factor 

^  c+(a) 

were  compared  with  those  of  the  authors  mentioned  above  for  the  same  values 
of  as  in  Figure  11.  It  is  obvious  that  there  is  an  almost  exact  corres¬ 

pondence  between  the  results.  The  above  authors  have  plotted  curves  for 
*  .9,  .8, . i,  and  although  it  is  not  shown  here,  the  equivalence 


displayed  in  Figure  14  remains  undiminished  for  each  and  every  value  of  the 


ratio  of  crust  thicknesses. 
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V. 

CONCLUSIONS 

The  conclusions  drawn  here  will  be  presented  in  the  order 
of  decreasing  generality. 

Although  the  first  of  these  conclusions  cannot  be  drawn  from 
the  material  in  this  thesis  alone,  it  is  included  as  an  indication  of 
the  efficacy  of  the  Wiener-Hopf  technique. 

In  the  preliminary  research  connected  with  this  work,  exposure 
to  literature  including  references  1,  4,  6,  10  and  11  has  shown  that  the 
type  of  problem  amenable  by  the  Wiener-Hopf  Technique  often  contains  com¬ 
plicated  mixed  boundary  conditions.  The  range  of  this  complicated  problem 
is  a  significant  extension  beyond  that  of  the  simpler  problem  for  which 
solutions  can  be  obtained  by  the  method  of  separate  e  v^iables.  Appli¬ 
cation  of  the  Wiener-Hopf  theory  is  facilated  by  variations  within  the 
technique  and  modifications  such  as  the  one  shown  in  Figure  2.  A  typical 
mixed  boundary  value  problem  is  that  of  satisfying  the  steady-state  wave 
equation  in  free  space  when  semi-infinite  boundaries  are  present. 

The  geometry  of  the  problem  considered  in  this  thesis  is  such 
that  the  Wiener-Hopf  technique  is  most  appropriate  for  obtaining  the  trans¬ 
mission  and  reflection  coefficients.  Aside  from  the  limitation  of  absolute 
values  imposed  on  the  coefficients  ar.J  the  ostensibly  negligible  energy 
lost  via'.;  propagation  down  the  duct  of  Figure  2,  the  analysis  produces 
accurate  results  for  a  large  range  of  crust  thicknesses  and  discontinuity 
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sizes.  A  solution  by  Sato  (reference  11)  and  a  modified  Green's  Theorem 
approach  by  Knopoff  and  Hudson  (reference  7)  both  entail  assumptions  which 
either  limit  the  range  of  results  or  make  them  approximate. 

A  key  to  evaluation  of  the  results  of  this  thesis  l'es  in  the 
comparisons  of  Figures  11  and  14.  Deviations  from  the  dashed  curves  of 
Figure  11  are  not  important  in  themselves,  for  as  is  shown  in  the  discus¬ 
sion  immediately  following  the  introduction  of  this  figure,  these  curves 
are  approximate.  Furthermore,  it  is  suspected  that  the  results  of  Knopoff 
ana  Hudson  are  too  large  for  intermediate  crust  thicknesses.  As  a  check, 
a  calculation  was  performed  of  TFT/TF1  using  the  amplitude  coefficients 

jT^j  and  Jt^J  corresponding  to  Knopoff's  and  Hudson's  approximation  kith 
+ 

^ assumed  to  be  unity).  Since  the  form  of  the  transmitted  Love  wave 
o  (ot) 

in  the  work  of  Knopoff  and  Hudson  is  equivalent  to  (72)-(74),  this  calcu¬ 
lation  is  a  test  of  the  accuracy  of  their  results.  The  resulting  values 
of  the  energy  transmission  coefficient  TFT/TFI  were  larger  than  unity  for 

.05  <  k  (h+H)  <  .38,  thus  labeling  Knopoff's  and  Hudson's  amplitude 
P1 

transmission  coefficient  as  t.oo  large  for  at  least  this  small  range  of 
crust  thickness. 

The  approximate  nature  of  Knopoff's  and  Hudson's  results  creates 

the  possibility  that  they  are  too  large  for  k  (h+H)  >  .38.  Strengthening 

31 

this  possibility  is  the  fact  that  oversize  coefficients  may  not  violate 

energy  restrictions  for  increasing  intermediate  crust  thickness  due  to  the 

decreasing  behavior  of  TFT/TFI.  Therefore,  the  fact  that  the  coefficients 

of  Knopoff  and  Hudson  are  larger  than  those  of  this  analysis  for  k  (h+H)<  6 

31 

must  in  the  main  be  attributed  to  error  on  their  part. 
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The  solution  to  the  problem  considered  in  this  thesis  is  in¬ 
complete  insofar  as  the  SH  body  waves  associated  with  the  branch  cut 
integral  in  (69)  have  not  been  discussed.  A  knowledge  of  the  energy 
contained  in  this  scattered  field  would  provide  an  absolute  check  on  the 
graphical  results  presented.  Although  the  comparison  with  Knopoff  and 
Hudson  reveals  the  accuracy  of  these  results,  knowledge  of  this  radiated 
energy  completes  the  picture  of  energy  apportionment.  This  picture  would 
be  a  criterion  for  determining  to  what  degree  the  modified  problem  of 
Figure  2  approximates  the  original  problem,  and  would  establish  the 
maximum  H  that  could  be  considered. 
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APPENDIX  A; 

FACTORIZATION  OF  £(v) 

Introductory  Remarks 

The  integral  formulas  (66) -(68)  are  as  yet  unable  to  produce  explicit 
expressions  for  the  total  field  for  they  contain  the  arbitrary  factors 
of  E(v).  There  is  a  general  method  of  factorization  (reference  11)  which 
allows  a  determination  of  e+(v)  and  a”(v)  at  least  as  a  product  of  ex¬ 
ponentials  (with  arguments  in  integral  form)  and  algebraic  terms.  If 
the  absolute  value  of  these  factors  is  all  that  need  to  be  determined, 
then  the  expressions  representing  them  can  be  reduced  to  a  form  suitable 
for  numerical  evaluation. 

The  method  is  based  on  the  theorem  that  a  function  q(v)  analytic 
in  the  strip  jlmuj  <  b  is  the  sum  of  two  functions,  one  being  analytic 
for  Imu  >  -  b,  the  other  for  Imu  <  b,  and  the  procedure  adopted  is  the 
application  of  Cauchy’s  Integral  Theorem  to  the  boundary  of  the  strip. 

The  integral  breaks  up  into  two  integrals  representing  the  functions 
mentioned  above.  In  general,  the  representation  is  as  follows: 


(A-l) 

where 


q(v)  ■  q  (v)  -  q  (v) 


|t(u)  ■  2^7  f 


ia  q(p) 


-®-  i  a 


(A-2) 


‘  1ST 


jo+ia 


0  -  17 


q(p) 


dp 


fc+i 


ia  p  -  u 


dp 
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q  (v)  is  analytic  for  Imv>-a,  while  q  (u)  is  analytic  for  Imv<a,  and 
thus  q(v)  is  analytic  in  the  st'^p  |lim?j<  b  since  a  can  be  chosen  as 
near  b  as  desired, 

Application  to  E(v) 

Equations  (A-l)  and  (A- 2)  can  be  employed  to  represent  Z(v) 


aS  rS  grated  in  (61) „  Sine  the  method  produces  a  sum  rather  than 

cr~(v)  + 

a  product,  this  is  equivalent  to  breaking  up  In  £(p)  into  In  a  (v)-  In  a  (v). 

This  may  be  done  as  long  as  In  £(u)  has  no  singularities  in  the  strip. 

Since  E(kj,  is  analytic  in  this  region  by  virtue  of  restriction  to  the 

]~2  ~ 

Reimann  s^eetwhere  R?  jv  *  k  >  0,  la  £(v)  possesses  the  same  property 

°2 

there  provided  that  £(v)  has  no  zeros  for  jlmi>|<|k  jsin  e.  Manipula- 

02 

tion  of  (60)  results  in  a  form  which  facilitates  the  location  of  the  roots 
of  E(v) .  This  form  is  as  follows. 


(A-3) 


E(v) 


1-  tan  s^H  tan  s^h 

r^2s2 

s.  tan  s.ri! - -  tan  s 

I  l  i  S 

11 


T tan  [s  (h  +  H)]  1 

1  U.s,  J 


r’l 


»Vl 


Roots  1  -  tan  s^H  tan  s^h  =  0  imply  that 


(A-4)  s^h  +  K)  *  n/2,3n/2, 


which  is  incons is tent  with  the  restriction  that  the  nature  of  the  period 

equation  places  upon  this  parameter.  From  Figure  3  it  is  e’  ident  that 

for  a  finite  layer  thickness  s^(h  +  H)  must  be  less  than  ^  ,  n  -  1,2  ... 

for  v  m  v  (oeriod  equation  root),  n  *  1,2,..,,  respectively.  Roots  K»~ 
IT  c 

of  the  period  equation 
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(A-5)  tan  !~s^(h  +  H)"|  - 


^2S2 

^1S1 


for  the  thicker  layer  have  been  assumed  to  lie  on  the  edge  of  the  ana¬ 
lytic  s~rip.  Therefore  E(v)  has  no  zeros  in  the  strip. 

In  order  to  insure  the  convergence  of  the  integrals  (A-2) , 

fi?  must  be  zero;  For  the  case  being  considered,  this  requirement 

IP!  "*  ® 

takes  the  form  £im  E(v)  =  +1.  From  (60),  with  s^,  s^  behaving  as  iv,  v 

M  -  » 


respectively,  as  Jvj  —  co,  it  can  be  shown  that 


(A-6) 


E(v) 


jvj-*  «  j  vj-oo 


ns  long  as  s^H  is  not  equal  to  zero. 

By  virtue  of  (A-6),  the  function  q  cannot  be  the  logarithm  of  E(v) 
alone,  but  factors  can  be  introduced  to  develop  a  function  that  does  satisfy 
the  necessary  limit  condition.  Consider  the  function 


(A-7)  P(v)  *  r  /v2  +  |k  j2sin2e  E(v). 

1  h 


1/2  22 

Obviously,  the  term  -  /  v  +  jk  J  sin  e  serves  a  two-fold  purpose;  it 
cancels  the  limitirg  1/v  dependence  of  E(v)  without  introducing  singularities 
in  the  strip  |lmv|<|k  j  sin  e  and  brings  the  magnitude  of  P(v)  to  unity 
a»  jyj-  oo.  Therefore,  the  method  can  be  applied  to  #n  P(v)  as  follows: 


(A-8)  q(v)  **  £n  °(v)  -  q"(v)  -  q  (v) 


*  -  ia 


(A-9) 


q+(v) 


llL dy 
P-V 


0  <  a  <  Jk  js in  e 
B2 


-  -  ia 
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(A- 10) 

oo  +  ia 

1  r  In  P  (0) 

q  U)  2^1  +  ,  T  -  »  ' 

-06  +  ia 

so  that 

(A- 11) 

-  eq‘CB)"  -  iv/7T 

or 

(A- 12) 

T.M - 2- - eq'C,,) 

/  2  |,  ,2  2 
\'V  +  ‘  32‘  Sin  G 

By  (61), 

(A-9)  ,  (A- 10),  and  (A- 12), 

(A- 13) 

+  ,  .  2 

a  (v)  *  exp 

Jv  +  iIkB2^Sin  G 

(A- 14) 

ct  (v)  =  Jy-i  |k  jsin  e  exp 

▼  On 

0  <  a  <  kfi  sin  e 


lUtf 


-ao  -  ia 


ia  -^n  P(p) 

-  dp  r 

p-y  M  J 


r<x>  +  ia  in  P(p)  ^ 

-  dp  1 

3-u  J 


-ao  +  ia 


The  functions  a  (v)  and  a  (v)  are  regular  and  teroless  in  the 
half  planes  Imy  >  -  a  and  Imv  <  a  respectively,  and  it  can  be  shown  that 
|cr  (v)  j  and  |ct  (v)  |  lie  between  positive  bounds  as  |v|  becomes  large  in 
their  respective  half  planes  of  regularity.  Thus  a  (v)  and  a  (v)  satisfy 
the  propert  es  listed  it  (61). 

A  calculation  of  the  exponentials  in  (A- 13)  and  (A- 14)  is  un~ 
likely  due  to  the  complexity  of  the  integrals.  In  order  to  benefit  from 
the  factorization,  the  integrals  in  (A- 13)  and  (A- 14)  must  be  reduced  to 
forms  that  are  solvable,  either  analytically  or  numerically. 

A  formal  manipulation  of  a  (u)  including  a  substitution  p™-o' 
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and  a  reversal  of  limits  of  integration  yields: 


(A- 15)  a  (-v)  -2  J 


■v-i|k  jsin  e 

°2  _  #  1 

v+i I k  i iin  e  +,  N 
1  82  a  (v) 


In  the  expressions  (75),  (76),  and  (79),  all  terms  are  real 
/-a-i|k  |  sin  c 

+  +  /  “  2 

except  cr  (K),  a  (a)*  and  J  —  sin  e  •  Therefore,  to  obtain  the 

V  a  ®2 

magnitude  of  the  transmission  and  reflection  coefficients,  one  needs  to 
take  the  absolute  value  of  only  these  factors,  all  others  presenting  no 
difficulty.  It  is  obvious  that  the  magnitude  of  the  above  radical  is 
unity.  As  for  |a  (K) | ,  the  analysis  is  as  follows:  By  (A-13), 


(A- 16)  |ct+(K) |  - 


/ K  +  ijk  Jsin  e 
”2 


-oo  -  ia 


U  +  i|ko  lsin  e 


{■-  w  r  '■  w  -I  )• 


•oo  -  ia 


Letting  the  width  of  the  analytic  strip  approach  zero  (e-*  0)  such  that 


(A-17)  J  K  +  i  Jk  jsin  e  /  K  ,  P(v)  -»  y  v  E(v) 


and  the  contour  of  integration  becomes  the  real  axis  with  suitable  indents- 
tions  around  P  »  ±  K,  ±  a,  and  ±  k  ,  and  noting  that  P(p)  is  even  yields: 


(A- 18)  |c+(K)|  -  ~z  exp  {  Im  [  "7  f 

V/K  L  L  J0 


p  E(P) ] 


0  p2  -  K2 
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An  inspection  of  (60)  shows  that,  for  0  >  V.  ,  p  real,  £(P) 

S2 

is  real,  so  that  the  integration  in  the  region  kQ  <  p  <  <»  produces  no 

w7 

imaginary  part.  For  0  £  P  <  k  ,  £(p)  is  eomplex^with 

02 


M  S 

u  “  — ~  sin  s.h  1m  [6(P>]-  cos  s.h  Re[6(p)] 
cot  s.H  pi.  s.  1  1 

(A- 19)  Re[£(o)]«- - —  +  — — - 


|Re[6(c)]|2  +  -  Im[6(p)]j^ 


and 


U-20)  Im[S(P)3  - 


s2 


-  sin  s.h  Re[6(p)]  +  cos  s  h  Im[$(p)J 

^1  1  _ l _ 

{Re[6(p)]}2  +  {ltn[6(P)]}2 


where 


^9  - 


ReCfiCp)]  “  s^in  Sjh  ,  IrnUCp)]  =  ~  s2  cos  si 


(A-21) 


;x  -  s^p)  »  Jk 


2  2 

Bf  P  ’  S2  ~  S2 


(p)  =  -  i  Jk 


2 - Z  - 

B2  '  P  -iS2 


Expressing  in[y  pl(p)]  as  4n[|~  pl(p)|]  +  i<Kp)  , where 


(A- 22)  $ 


--•MSW} 


(A- 18)  can  be  written 


(A-23)  |/(K)|  --|exp{.Sj  02  ^Mdp} 

a/*  0  P  -K 


lit)'!- 
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This  form  of  the  factor  is  still  generally  unsolvable,  for 
the  values  of  k  and  K  are  usually  unknowns,  and,  even  if  they  were 

b2 

available,  the  complexity  of  the  integrand  makes  any  attempt  at  formal 
integration  futile. 

It  is  possible  to  make  a  change  of  variable  that  surmounts 

these  difficulties.  Let  p  *  k  T)  which  results  in 

P2 


(A- 24) 


where 


2 


dT) 

2 


(A-25) 


Sl(7}) 


The  ratios  of  propagation  constants  can  be  obtained  from  the  period 
equation  and  knowledge  of  the  value  of  the  ratio  of  the  shear  wave 
velocities  3^  and  32,  making  calculation  of  the  real  and  imaginary 
parts  of  E(T1)  possible.  Although  analytical  integration  is  still 
fairly  impossible,  (A-24)  is  of  a  form  such  that  |ct  (K) ]  can  be 
evaluated  numerically. 

By  a  similar  procedure  the  factor  |a  (a)|  is  represented  by 
expression  identical  to  (A-24),  except  that  a  replaces  K.  Thus  the 
expressions  (75), (76),  and  (79)  for  the  transmission  and  reflection 
coefficients  can  be  programmed  on  a  digital  computer,  the  integration 
in  (A-24)  being  approximated  numerically.  The  end  results  of  these 
calculations  are  in  principle  the  coefficients  desired. 
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FIGURE  4  FUNDAMENTAL  LOVE  WAVE  PHASE  VELOCITY 


branch  cut 


-  PLANE  CONTOUR  FOR  TRANSMITTED  FIELDS 
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FIGURE  8  THE  TRANSMISSION  COEFFICIENT  |T2| 
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FIGURE  9  THE  REFLECTION  COEFFICIENT  Rz 


IN  TKE  TRANSMITTED  LOVE  WAVE 


FIGURE  II  COMPARISON  OF  |TT|  WITH  THE 

RESULTS  OF  KNOPOFF  AND  HUDSON 
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FIGURE  12  THE  COEFFICIENT  |T,|  VIA  THE 

METHOD  OF  KNOPOFF  AND  HUDSON 
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FIGURE  14  COMPARISON  OF  APPROXIMATE  |TT|  WITH 
THE  RESULTS  OF  KNOPOFF  AND  HUDSON 


